CLNS 01/1750 



Integrable Boundary Conditions and Reflection Matrices for the 

0{N) Nonlinear Sigma Model 



M. Moriconi ^ 

Newman Laboratory of Nuclear Studies, Cornell University 
Ithaca, New York 14853, USA 



Abstract 

We find new integrable boundary conditions, depending on a free parameter g, 
for the 0{N) nonlinear a model, which are of nondiagonal type, that is, particles 
can change their "flavor" through scattering off the boundary. These boundary 
conditions are derived from a microscopic boundary lagrangian, which is used to 
establish their integrability, and exhibit integrable flows between diagonal boundary 
conditions investigated previously. We solve the boundary Yang-Baxter equation, 
connect these solutions to the boundary conditions, and examine the corresponding 
integrable flows. 
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1 Introduction 

One of the main problems in boundary intcgrablc field theories is to find boundary con- 
ditions that preserve the integr ability of a given bulk integrable field theory, that is, 
integrable boundary conditions. This has been carried out for several models by direct 
inspection of higher-spin conserved charges [1, 2, 3, 4]. In [3] (see also [5, 6]) it has been 
shown that models with 0{N) or SO{N) symmetry admit a family of A'^ -|- 1 integrable 
boundary conditions, which are obtained by requiring n field components to satisfy Neu- 
mann and the remaining N — n to satisfy Dirichlet, n — 0,1, . . . , N. These conditions 
are of diagonal type, which means that scattering off the boundary will not change the 
fiavor of the incoming particle. We will refer to this type of boundary scattering prob- 
lem as the diagonal case. A natural question one would ask then is if there are more 
integrable boundary conditions for these models with 0{N) symmetry, or if the diagonal 
case exhaust all possibilities. 

The purpose of this paper is to study the possible integrable boundary conditions for 
models with 0{N) global symmetry. We find nondiagonal integrable boundary condi- 
tions for the 0{N) nonhnear a (nlu) model. These are genuinely nondiagonal boundary 
conditions in the sense that they can not be put into diagonal form by a redefinition 
(global rotation) of the fields. These conclusions should carry out to other models, like 
the SO{N) Gross-Neveu and principal chiral models. 

The structure of this paper is as follows. In the next section we review some of the 
generalities of boundary integrable field theories which are relevant for our discussion. We 
also review the Goldschmidt-Witten argument, and use it to establish the integrability of 
the 0{N) nonlinear a model. In section 3 we look at a very simple example, to illustrate 
some of the issues that are relevant in our analysis. In section 4 we introduce the new 
integrable boundary conditions for the nonlinear a model, which are the main concern 
of this paper. In section 5 we compute the nondiagonal reflection matrices associated to 
these boundary conditions by solving the boundary Yang-Baxter equation, in particular 
we look at the 0(2) case, and present a solution for the 0{2N) case to which we do not 
know what are the corresponding boundary conditions. Finally in section 6 we present 
our conclusions and some possible directions for further work. 

2 Integrable Boundary Field Theory 

Let us consider a bulk integrable field theory defined by an action 




(2.1) 
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where CB{4>,d(j)) is the bulk Lagrangian, which depends locally on the field and its 
partial derivatives ^. We say that it may be defined by a local action because there are 
several bulk integrable models that are not, such as perturbed conformal field theories. 

As is well known, in two-dimensions integrability has several important consequences: 
in a multi-particle scattering process there is no particle production and the set on in 
and out momenta is the same, the S'-matrix factorizes into a product of two-body S*- 
matrices, and the two-body S'-matrix satisfies the Yang-Baxter equation, besides the 
usual requirements of unitarity and crossing-symmetry. These conditions, together with 
the knowledge of the symmetries of the S'-matrix, often allow one to determine the two- 
body S'-matrix exactly, up to CDD factors. 

When we consider this theory on the half-line, the introduction of a boundary may 
easily break the conservation of bulk charges, as it is seen in the simple case of linear 
momentum. Therefore we are not assured that the resulting model will be integrable, 
and we must specify boundary conditions. These boundary conditions correspond to 
boundary dynamics, which may be described by some local lagrangian. The boundary 
version of a given bulk model is defined by the action 

/O r+oo r+oo 

dxi / dxo jCb{(P, d(t>) + / dx^ A(0L=o ' _n) ' (2-2) 
-oo J —oo J —oo — ^ 

where Cb{(f>\^=Q ,0 ^) is the boundary action, depending only on the values of the field 
(f) and its time derivatives, at a; = ^. The boundary Lagrangian defines the boundary 
dynamics, and it's easy to see that it induces boundary conditions through equations of 
motion. For example, for a scalar boson we have 



Ol<P\x=0 



d(j) 



+ Oq 



x=0 



d(f) 



= . (2.3) 

x=0 



In order to have a boundary integrable model, we have to pick a suitable boundary 
action such that a combination of bulk conserved charges is still conserved. In order to 
discuss boundary conserved charges, we briefly review the bulk case. 

A local (bulk) conservation law of spin s is written as 

9_ J^+') = d+R'f-^'> and d+A'+'^ = , (2.4) 

and from these equations, one easily establishes the conservation of the following charges 
Q+ = dxi ( Ji'+'^ - R^r^^) and Q_ = dx^ (A'^'^ - R^^'^) (2.5) 

J —oo J —oo 



^Here the field (j) is symbolic and may denote a collection of fields, bosonic or fermionic fields and so 

on. 

^We use an upper dot for time derivative at the boundary. 
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In proving that these charges are conserved we have to use the fact that we can discard 
surface terms. When we restrict our model to the half-hne we can not do that with the 
surface term at x — 0. On the other hand, if the following condition [1] is satisfied 

ji_s+i) _ ^(s+1) ^ _ 1^^^ ^ d^^^^ ^2.6) 

CLv 

where E(i) is a local field, then 

Q = f dx, (Ji^+'^ + 4^'^ - R^r'^ - - (2.7) 

J—oo 

is a conserved charge. For example, in light-cone coordinates the energy-momentum 
charges are Q± — E±P, but after the introduction of a boundary only Q — Q++Q_ — 2E 
is conserved, with T,(t) — 0. 

We should look for the possible boundary conditions that we can impose to a given 
model such that 2.6 holds. Once this is done one should look for solutions of the boundary 
Yang-Baxter equation (bYBe) which can be related to the proposed boundary conditions. 
We describe the bYBe in the next subsection 

2.1 The Scattering Matrix: Generalities 

Boundary integrability is encoded in the boundary Yang-Baxter equation (bYBe) in a 
similar way as bulk integrabihty is encoded in the Yang-Baxter equation (YBe). We 
should stress that, unlike the bulk case, we do not know what is the symmetry group 
at the boundary a priori. In the bulk case the knowledge of the bulk global symmetries 
greatly simplify the form of the jS'-matrix. Unfortunately we do not have a similar set 
up for the boundary case, and therefore the analysis of the possible structure of the 
reflection matrix will have to rely on general physical arguments. Recently this issue has 
been addressed in [10] for the case of the nonlinear Schroedinger model. We review now 
a few well-known concepts. For a more complete discussion, see [13] 

The Hilbert space of in (out) asymptotic states is spanned by the multi-particle states 

\Ai^{ei),Ai,{e2) . . . A„(^n))m(out) = ^^(^1)^2(^2) • • • A„ (^n) |0)in(out) , (2.8) 

where the {Ai{9)} are the Faddeev-Zamolodchikov (FZ) operators that create the one- 
particle asymptotic states, the {9i} are the rapidities and > > . . . > for 
in-states and the other way around for out-states, and \0)in{out) is the in (out) vacuum. 
We assume henceforth that the in and out vacuum are the same. 

Multiparticle scattering processes factorize in a product of two-body yS'-matrix, which 
is defined in terms of the FZ operators by 

Ai{9i)A,{92) = S^;{9, - 92)Ai{92)Ak{9i) . (2.9) 
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The YBe is obtained by requiring the associativity of the algebra defined by 2.9 

In defining the refiection matrix we have to change the Hilbert space. Boundary 

scattering processes are defined in terms of the boundary state |0)b — B\0), where B is 

the so-called boundary-state operator [1]. 

The Hilbert space is spanned by . . . Ai^{6n))B = ^n(^i) • • • b, where 

the in (out) states are obtained by setting 9i > ... > 6n > {6i < ... < On < 0). 

Similarly to the bulk case, the refiection matrix is defined by 

Ai{e)B = R{{e)Aj{-e)B . (2.10) 

We represent the ^-matrix and the reflection matrix graphically as in figure 1. 



1, ,k 




i j 



Fig. 1 The S-matrix and the refiection matrix 
The bYBc is obtained by requiring the compatibility of 2.10 with 2.9 by looking at 




Fig. 2 The boundary Yang-Baxter equation. 
In terms of the two-body ^-matrix and refiection matrix, it reads 

Sr{0)K{0i)SU(^+)R',{92) = Rm)S-{0+)RUOi)S'U(^) , (2.11) 
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where we introduced the variables 6+ — 6i + 62 and 6 — 61—62, and sum over m,n,p 
and q. 

Besides the bYBe, one must impose unitary and crossing-symmetry conditions on the 
reflection matrix. We quote them here, and refer the reader to [1] for more details. The 
unitarity condition is 

i^iO)Ei{-6) = Si , (2.12) 
and the boundary crossing symmetry is 

RKj-0) = Si^i{26)Ri{^ + 6), (2.13) 

where we are assuming implicitly that the particles are invariant under charge conjuga- 
tion. For the general case we refer to [1] and the appendix of [12] . 

2.2 The 0{N) Nonlinear Sigma Model 

In this subsection we briefly review the main features of the 0{N) nonlinear sigma model. 
The nlcr model is defined by the following Lagrangian 

Cnia = :^dn-dn (2.14) 

where n is an A'"-dimensional vector subject to the constraint n-n = 1. This constraint can 
be implemented by means of a Lagrange multiplier A in the action. It is straightforward 
to see that the equation of motion for n in light-cone variables {x± — xi) is 

d+d-ii = —n{d+n ■ d-u) . (2-15) 

This will be used in establishing the integrabihty of the boundary conditions of section 
4. 

The main features of the nlcr model is that it is classically conformal, asymptotically 
free, and displays dynamical mass-generation. 

2.3 The Goldschmidt-Witten Argument 

We can use the fact that the nla model is classically conformally invariant in order to show 
that it is indeed an integrable model at the quantum level. Starting with the classical 
conformal symmetry and exploring the structure of the possible anomalies that appear 
after quantization, Goldschmidt and Witten [9], and earlier Polyakov [7], have shown that 
some conserved currents arc not spoiled by quantization, that is, there are anomalies, but 
their structure is such that one can rewrite the quantum corrections to the conservation 
law as a total derivative, allowing us to write down non-trivial conserved charges. This 
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argument can be used to establish the integrabihty of the Gross-Neveu model [8] and 
of the principal chiral model [9]. In this subsection we review the Goldschimdt-Witten 
argument and write down the conserved currents that will be of use to establish the 
integrabihty of the new boundary conditions for the nlo" model. 

The Goldschmidt-Wittcn (GW) argument starts by observing that the trace of the 

energy- momentum of a classically conformal-invariant theory vanishes, that is = 

TL-i- — 0. This implies that the conservation of energy-momentum reads 

= and = (2.16) 

which, on its turn has as a consequence that for any integer n 

= and a+(r__)" = o . (2.17) 

Therefore, we generate towers of classically conserved currents. 

After we quantize the theory the equations 2.17 do not make sense any more, since 
we are considering products of operators at the same point, and these operators should 
be redefined. In general these equations will be spoiled by anomalies. The right hand 
side of 2.17 will be nonzero, but there are still several constraints we should impose on 
the possible terms that are generated: whatever comes due to anomalies has to have the 
correct dimension, Lorentz weight ^ , and group theoretic properties as the left-hand side. 
By analyzing all the possible terms that may appear on the right-hand side of 2.17, with 
n — 2, GW showed that they can all be written as total derivatives, which implies that 
we have a nontrivial conserved current. 

For the nlcr model they showed that the possible anomalies can be written as 

a+(T__)2 = ci a+(a!n • a!n) + C2 a_(a+n • 9_n d.n ■ d.n) + cg d.{dtn ■ d+n) , (2.18) 

with some constants {q} and an analogous equation for d^{T )^. 

In [3] wc used these charges to establish the integrabihty of the aforementioned diag- 
onal boundary conditions. We will sec that we can use them to establish the integrabihty 
of nondiagonal boundary conditions too. 

2.4 The Exact 5-matrix 

Once the nlcr model has been established to be integrable we can proceed and compute 
its two-body ^'-matrix. This has been done in [11], and we quote it here for further 
reference. Since the bulk is 0{N) invariant, the (S-matrix has to be of the form 

Sf^{9) = ar^{9 )6^,6'' + a2{e)5^5\ + a^{e)^,5'; (2.19) 

■^If under a Lorentz transformation parametrized by rapidity a, a quantity 4> €^"0, we say that (j) 
has Lorentz weight p. 
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where the (Ti{6) are determined by the solving the YBe [11], unitarity, and crossing- 
symmetry. In [11] a large- expansion check has been performed to show that this is 
indeed the exact ^'-matrix of the 0{N) nla model. The ai{9) are given by 

<yi{0) = --^a^iO) , as{e) = -\2{e) and 

(Q\ _ V2 ^ 27r ^ 27r^ V ^ 2-k ! \2 2n ^ V 27r 27r/ /o oriN 

r(i + A-f)r(-f)r(i + i|)r(i + A+«)- 

where A = 27t/{N — 2). We will also need the isoscalar ^'-matrix element cri{9) — 
Nai{9) + (72 (^) -I- (73 (^), which can be written as 

C7i{e) _ g^g . (2.21) 

This will be useful when we solve the boundary crossing unitarity condition in section 5. 
For N — 2 one has to be careful, since A — > oo but (72 (^) 0. In this limit we get 

977]- /*77r 

CTi{0) = -. z/W ' ^2^ = , and as{e) = -^f{e) (2.22) 

where the function f[6) is given by 

^ ^ 27r/ V2 27r/ 

The 0{2) nl(7 model is not a simple massless free boson, as the map rii — cos(y^0) 
and n2 — sin(^0) would suggest. It is actually the sine-Gordon model at = Stt, 
the sine-Gordon potential being a marginally-relevant perturbation at this point, and it 
describes the Kosterlitz-Thouless point of the classical XY model. 



3 A Very Simple Example 

Before we analyze the nondiagonal boundary conditions for the nla model, let us look at 
a very simple example of nondiagonal scattering, consisting of 2 free bosons (pi and 02, 
with equal mass m, on the half-line, which captures some of the main physical aspects 
for the case of the nla model. The action we consider is 

/oo rO 1 1 1 ]^ roo 

dxo dxi -{d(f)iY + -mj(f)j + -{d(f)2T + -mlcpl + A(0i, 0i, 02, ^2) , 
-00 J -—00 ^ Z Z Z J —00 

(3.1) 

where is the boundary action, to be chosen shortly. 
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Clearly any quadratic form at the boundary can be solved exactly, being obviously 
integrable. We have consider two possibilitites ^ 



- 

'~'h - 



W2 



(3.2) 
(3.3) 



Recall that a boundary Lagrangian Cb{<Pi) induces the following boundary condition 

dC 



dx 



x=Q 



d(j)i 



(3.4) 



x=0 



Since the masses of the two particles are the same, which is required to have nondi- 
agonal scattering, we can define new fields r]± by 

1 



(01 ± 02) , 



(3.5) 



which leave the bulk action with the same form but takes the boundary action to 

Cl^g,{vl-rl) (3.6) 

and therefore we conclude that the boundary condition induced by 3.2 is not really 
nondiagonal. 

On the other hand, for the second boundary Lagrangian it is not possible to find an 
orthogonal transformation which diagonalizes the boundary Lagrangian. In this sense we 
can say that the boundary Lagrangian II corresponds to nondiagonal boundary scattering. 

It is fairly easy to find the reflection amplitudes for 3.3. The boundary conditions 
induced by 3.3 are (using g for 

d(f)2 



dx 
d(f)2 



= g 



dx 



x=0 



x=0 



dt 



dt 



x=0 



x=0 



By using the mode expansion for the fields, 

^.(^j. _ J ^Q^j^^^Q-^^{imcosh{0)t-imsmh{0)x) _|_ j'^jg(-imcosh(e»)t+imsinh(6')a;) j 



(3.7) 
(3.8) 

(3.9) 



where the Ai{9) are the FZ operators, and the definition of the reflection matrix Rl{9) 
given in 2.10, we find the scattering amplitudes to be 



( smh^(e)-g^cosh-^(6') gsinh(2f>) \ 

sinh^(6»)+ff2cosh^(6») {B)-\-sm!ci^ {&) 

-gsinli(26') sinlig (6')-g^ cosh^(6>) 



(3.10) 



"'The choice £(, = + ^202 is trivial, since there is no couphng between the two fields. We are also 
disregarding possible boundary mass terms. 
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Notice that, as expected, the diagonal amphtudes are the same, and the off-diagonal 
ones have opposite signs. Moreover if we expand the elements of the reflection matrix in 
powers of the coupling constant the diagonal elements are even functions of g, whereas 
the off-diagonal ones are odd functions of g. This can be understood by noticing that 
each flavor change at the boundary is accompanied by a factor of and so in order 
to have diagonal (nondiagonal) boundary scattering we need an even (odd) number of 
(yf's. We will use some of the intuition from this very simple model in the study of the 
nondiagonal boundary conditions for the nlcr model. 



4 Boundary Conditions for the nlcr Model 

As we have seen in section 3 there are essentially two possibilities we could try for the 
boundary conditions in the nlcr model. One of them, 3.2, is not even integrable in this 
case. Wc are left then with the boundary conditions induced by a boundary lagrangian 
of the form 3.3, 

Snia = Sb+ dxo-Mijriirij , (4.11) 

J —oo Z 

where we have to have Mij — —Mji, and the indices i and j run through a subset 
of {1,2,..., N}. For the remaining indices, which correspond to particles scattering 
diagonally, we choose Dirichlet boundary conditions 

donj\^=o = 0- (4.12) 

The reson for this choice will become clear in the next section. Prom 4.11 we get the 
following boundary conditions for the non-diagonally scattering fields 

^1 rii 1 00=0 = Mij do nj\oc=o ■ (4.1 3) 

We have to check that the condition 2.6 is satisfied for these boundary conditions, that 
is 

(9_n • 9_n)2 - (9+n • d+nf + ci(9^n • a^n - din ■ din) + 
-|-C2(9+n • 9_n 9+n • d+n — 9+n • 9_n 9_n • 9_n) -|- 

+C3{dln ■ d_n - din ■ d+n)\^_^ = ^E{t) (4.14) 

for some local field S(t). We discuss this condition term by term. The first term can be 
simplified to 

(9_n • (?_n)^ — (9+n • 9+n)^ = 8 no • ni (no • no + ni • ni) , (4-15) 
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where the subscripts and 1 denote time and space derivatives, respectively. The "ci" 
term can be simphfied to 

• d^n — d^_n. ■ d^n — 16 noo ■ noi — 8 no ■ ni (no ■ no — ni ■ ni) , (4.16) 

the "C2" term to 

d+n ■ d-U d+n ■ d+n — 9+n • 9_n 9_n • 9_n = 4 no • ni (no • no — ni ■ ni) (4-17) 
and the "03" term to 

d^ii ■ d-H — d^ii ■ 9+n = 8 9o(noi • no — noo • ni) — 2 no ■ ni (no ■ no — ni ■ ni) (4.18) 

We see from these that we have to have no • ni = in order that 2.6 is satisfied. Moreover, 
if we consider 4.12 and 4.13, we immediately see that the one "bad" remaining term 
16 noo • noi vanishes. This establishes the integrability of 4.12 and 4.13. 

Since M is an antisymmetric matrix, there is a real-orthogonal matrix O, such that 
O^MO is block-diagonal, each block being an antisymmetric 2x2 matrix. If M has k 
zero eigenvalues we readily see that its structure will have 0(2) x 0(2) x . . . x 0(2) x 0{k) 
symmetry, where there are {N — k)/2 0(2)'s (and therefore N — k must be even). The 
0(2)-symmetric blocks correspond to interactions between pairs of field components at 
the boundary, that is, a particle of a given type, say i, can scatter diagonally or onto 
some other fixed particle of type j, where the pairing is fixed. 

The question we should try to answer now is, in which ways can we break the boundary 
symmetry? In the diagonal case, the bulk 0{N) symmetry is broken to 0{k) x 0{N — k) 
at the boundary, with k = 0,1, . . . , N. We will be looking at solutions that respect this 
symmetry, and so we should consider matrices M symmetric under 0{k) x 0{N — k). 
But this implies that the only case we should consider is 0(2) x 0{N — 2), where the 
0(2)-symmetric part corresponds to non-diagonal scattering, and the remaining 0{N — 2) 
to diagonal scattering, that is, the components rii satisfy 

dini\a:=o = gdon2\x=o , din2\x=o = -gdoni\x=o 

9on, U=o = for j = 3, 4, . . . , AT . (4.19) 

We have chosen the diagonally-scattering field components to satisfy Dirichlct boundary 
conditions, but at this point it seems that wc could have chosen them to satisfy Neumann. 
As we will see later, the solutions of the boundary Yang-Baxter equations are such that 
in the "diagonal" limits the correct choice for the diagonally-scattering components is 
indeed Dirichlet boundary conditions. 

In [14] Corrigan and Sheng have established the classical integrability of the Neumann 
boundary condition for the 0{N) nla model, for all N, and found one extra boundary 
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condition for N = 3. In [3] we showed that this special boundary condition was integrable 
at the quantum level to. We will show now, that it is actually related to the ones we 
propose in this paper. 

The explicit form of the Corrigan-Sheng boundary condition for the 0(3) nlcr model 

is 

9in|j;=o = X 9on|a;=o + (n • /c X 9on) n|j;=o and k ■ don\x=o ^ (4.20) 

where k is an arbitrary vector. Let us choose k = g (0,0,1), where g is a coupling 
constant. This means that don^ = and so is a constant at the boundary. This, 
together with the constraint n\+ n\ + n\ = 1, imply n\ + n\ = constant at x = 0. The 
Corrigan-Sheng boundary condition reduces to 

9ini|j;=o = g dQn2\x=o , din2\x=o ^ -g doni\a:=o 

dons\x=o = . (4.21) 

This is the same boundary condition we are considering in the 0(3) nlcr model case. 

5 Reflection Matrices 

In this section we will compute the reflection matrices associated to the boundary con- 
dition identified previously. We will concentrate on the nlc model, since the ones for 
the GN and PCM models, can be found by attaching suitable CDD factors. Since we 
are interested in the integrable flows for these boundary conditions, we will revisit the 
diagonal case before. 

5.1 The Diagonal Case 

In [15] the reflection matrices associated to the diagonal boundary conditions found in 
[3] were computed (see also [16]). One question one may ask is if there are integrable 
diagonal flows, that is, if it is possible to go from a conflguration with, say, M components 
satisfying Neumann and the remaining N — M Dirichlet boundary conditions, to some 
other other configuration with M' Neumann and N — M' Dirichlet boundary conditions. 
We will see that this is not possible. 

Let us suppose that there are two flavors i and j that scatter diagonally. Consider the 
scattering process \Ai{6i)Aj{62)) — > \Ai{62)Aj{6i)). Let's call the reflection amplitude for 
the i particle Ri{0), and for the j particle R2{0) The bYBe for this process is, therefore 

RM)R2{e2)a2ie+)asie_) + R2{02)R2{0i)asie+)a2ie_) = 

= Ri{ei)Ri{e2)a2{e.)a3{e+) + R,{e2)R2{e,)as{e.)a2{e+) . (5.22) 
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Dividing by i?2(^i)i?2(^2)(72(^+)(72(^-) and introducing X{e) = Ri{e)/R2{e), and Si{e) = 
(^i{0)/cr2{0), we get 



By taking the limit 02 ^ di, we arrive at the following differential equation 

X^{0) 



(5.23) 



1 



20 



(5.24) 



The solutions for 5.24 are X{0) — 1 which means that the amphtudes Ri — R2, or 



c + 



(5.25) 



Notice that this result is true for any pair of diagonally scattering flavors. Using 5.25 we 
will show now that there are no integrable diagonal flows. 

Suppose there are 3 diagonally scattering flavors, and call them 1, 2 and 3. From 5.25 
we know that 

Ri{e) c-0 Ri{e) c'-0 . R2{0) c"-e 



and 



(5.26) 



R2{0) c + ' Rs{0) d + R^{0) c" + 

with some constants c, c' and c". But these equations are incompatible, imless one of the 
ratios is 1. This proves that we can not have more than 2 different types of reflection 
amplitudes in the diagonal case. Moreover, as it was shown in [3], if we assume that 
there arc only 2 types of reflection amplitudes, the constant c in 5.24 can be fixed using 
the bYBe to be c = -i^^^- 

Now we can easily show that there are no diagonal integrable flows. Since the ratios 
are flxed by the bYBe, we simply can not go diagonally from a reflection matrix where 
there are k particles satisfying Neumann boundary condition, to one where there are 
k' ^ k. 



5.2 The Nondiagonal Case 

In this subsection we use the bYBe 2.11 in order to find the refiection matrices corre- 
sponding to the boundary conditions identified earlier. Using the intuition gained in 
section 3 we have the following ansatz for the refiection matrix 



/ A{0) B{0) 
-B{0) A{0) 



R 
















R^{0) 








R^{0) 



\ 



(5.27) 
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where A, B and Rq are to be determined. 

Let us consider the scattering process \Ai{6i)Ai{62)) — ^ ^2)), where i is 

any of the diagonally scattering particles, that is i > 2 The bYBe for this process is 

R^{e^)R^{e2)a^{e+)a2{e) + R^{e2)A{ei)a2{e+)a^{e) = Ro{e,)Ro{e2)as{e)a2{e+) + 
+A{e,)A{e2)a2{e)as{e+) - B{e,)B{e2)a2{e)as{e+) . (5.28) 

Dividing this equation by Ro{9i)Ro{92)a2{9+)a2{9), and defining X{9) = A{9)/Ro{9) and 
y(^) = S(^)/i?o(^),weget 

{X{9,) - X{92))d3{9) = {X{9i)X{92) - 1)S3{9+) - Y{9i)Y{92)S3{9+) , (5.29) 

and by taking the limit 9i — 6*2, we obtain the following differential equation 

Notice that if Y{9) = this equation reduces to the one for the diagonal case 5.24. Since 
we have two unknown functions now, we need one more equation. This is accomplished 
by the bYBe for |^i(^i)A(^2)) ^ |A(-^i)^(-^2)), which reads 

Ro{e2)B{e,)a2{e+)a3{e) = R,{dr)B{92)a2{9+)a^{9) + 
+A{9^)B{e2)a2{9)a^{9+) + B{9^)A{92)(J2{9)a^{9+) . (5.31) 

After dividing by i?o(^i)-Ro(^2)c2(^)c2(^+), and taking the hmit 9i — > ^2, we get 

^.y^e) ^ m . (5.32) 

We can easily solve equations 5.30 and 5.32 by introducing Z±{9) = X{9) ±iY{9), which 
satisfy 

d9^^^^^ - 29 ■ ^^-^^^ 
This is the same equation we obtained earher for the diagonal case 5.24. The solution 
for X{9) and Y{9) are, therefore 

i(c-e c'-9\ ^ 1 fc-e d-9\ 

Xi9) = - H and Yi9) = — (5.34) 

^' 2\c^9 c' + 9) 2i\c + 9 c' + 9j ^ ' 

Notice that the purely diagonal case corresponds to Y(6) = 0, that is c = c', and we 
recover 5.25. One aspect that may look puzzling is the fact that, at this point, we have 
two constants to be fixed, c and c', but only one coupling constant g. 
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In order to find a further constraint involving c and c', consider the scattering process 
|A(^i)^i(^2)) ^ |Ai(-^i)^(-^2)). The bYBe reads 

A{e^)B{e2)Me+)a:,{e) + cj2{e)a2{e+)) + A{e2)B{e^){ne+)a^{e) - a^{e+)a2{e)) = 

(AT - 2)R{e^)B{e2)ai{e)a^{e+)) (5.35) 

After dividing by R{9i)R{92)<J2{d)<^2{d+)-, we plug the solutions 5.34 into 5.35 to obtain, 
after an elementary but tedious calculation, that 

- 4 

c + c' = -i7r^^ (5.36) 

We have checked that the expressions 5.34 together with 5.36 solve all the remaining 
bYBe's. 

There are two ways to get diagonal reflection matrices from 5.34 and the constraint 
5.36. The first one is by setting c = c', and we obtain the result for the diagonal case 
with 2 fields satisfying Neumann boudary condition and the remaining N — 2 Dirichlet 
[3]. This is the deciding factor that made us choose Dirichlet boundary conditions for 
the diagonally-scattering fields in section 4. The other way to have diagonal refiection 
matrices is by having both c and d go to cxd. In this case we have X{9) = 1 and we 
conclude that all field components satisfy the same boundary condition, that is, Dirichlet. 

Therefore if we write c = -iTT{N - 4)/(2iV - 4) + ^{g) and c' = -iTr{N - A)/{2N - 
4) — ^{g)i where ^{g) is a (unknown) function of the boundary coupling constant, we 
have that, for = (2 Neumann, N — 2 Dirichlet) we should have {(0) = (and so 
X{e) = {-i7r{N - 4)/(2iV - 4) - e)/{-m{N - 4)/(2iV - 4) + 6)), and for ^ ^ oo (all 
Dirichlet), X{9) — 1 and so ^(oo) = oo. We see, then, that there are integrable flows 
between diagonal reflection matrices, and that the reflection matrices we found have 
the correct behaviour, as compared to the one dictated by the microscopic boundary 
Lagrangian. 

All that is left to do now is to solve the boundary unitarity and boundary crossing- 
symmetry conditions, which will be dealt with shortly. 

Based on the "very simple example" studied before and the solution we found for the 
case of 2 nondiagonal flelds, we are tempted at constructing other solutions of the bYBe 
in a similar fashion as the one just done. As we will see now, there are no such solutions. 
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Consider the following ansatz for the 0(5) model 

/ A{e) B{e) 

-B{e) A{e) 







m 



\ 




Ro{e) J 

This reflection matrix should correspond to the following boundary action 





-B{0) A{9) 




(5.37) 



(5.38) 



In principle we could have two different coupling constants, one for the fields ni and 712 
and another for and 714, but we are choosing them to be equal in order to simplify the 
argument. Since this case is a special point of the general case, the absence of solutions in 
this case will establish our claim that there are no solutions with more than 2 nondiagonal 
fields. We will look at the 0(5) case, but the argument is the same for 0{N) with more 
than on nondiagonal block. 

The bYBe for the processes |>li(^i)A(^2)) ^ |A(-^i)^2(-^2)) and 1^1(^1)^5(^2)) 
— > |^5(— ^i)^2(— ^2)) are essentially the same as in the case for 2 nondiagonal fields only, 
and so the ratios A/Rq and B/Rq have the same form as the ones given in 5.34. If we 
consider now the bYBe for the process 1745(^1)745(^2)) — ^ |^i(— ^i)^i(— ^2)) we find that 
there is no solution besides the diagonal one (c = c') for the constants c and c', and 
therefore there is no solution of the bYBe consistent with 5.37. 

Before we proceed, we should mention another solution of the bYBe for even N , but 
to which we do not know how to assign microscopic boundary conditions. 

Let us consider the following block-diagonal ansatz for the refiection matrix 



m B{0) 

-B{e) A{e) 










m B{e) 

-B{9) A{9) 



V 



(5.39) 



The bYBe for the process 1^1(^1)^3(^2)) ^ |A3(-^i)^2(-^2)) is 

A{92)B{9,)a2{9+)as{9) = A{9^)Bi92)a2i9+)a■,i9) + 

+ {A{9,)B{92)+A{92)B{9,)) a2{9)a^{9+) . (5.40) 

Dividing this equation by (72(^+)cr2(6')74(6'i)74(6'2), and taking the limit 9-i — > ^2, we get 

d (B{9)\ 1 (B{9)\ 



d9 \A{9) y 9 \A{9) 



(5.41) 
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whose solution is 



m 

A{9) 



— a 9 , 



(5.42) 



where a is a constant. We have checked that this alone solves all the remaining bYBe's, 
and therefore a is left as a free-parameter. All we need to do is to fix A{9) using unitarity 
and boundary crossing symmetry. Since their solution is quite straightforward, we will 
not quote it here. The only diagonal reflection matrix that can be obtained from this 
solution is by taking a — 0, which corresponds to all components satisfying Neumann 
boundary conditions. If a — > oo we have a block diagonal solution, each block being 
antisymmetric. Note that one must be careful in taking a — > oo, since unitarity and 
crossing symmetry will give a pref actor with an overall dependence on a, which precisely 
cancels the a dependence of B{9) in this hmit. 



5.3 The 0(2) Case 

In the limit ^ 2 we can solve the bYBe completely. In this subsection we will not 
write down the bYBe's explicitly, since they can be derived quite easily, and the 0(2) 
case is not our main interest. We will describe which equations one needs to solve to fix 
the structure of the refiection matrix completely, in any case. 
Initially consider the general form for the refiection matrix 



m 



f{9) w{9) 
h{9) g{9) 



The bYBe for | ^(^1)^1(^2)) ^ | A(-^i)^i(-^2)) gives 

K9r)w{92) = h{92)w{9,) . 
If these amplitudes are non-zero, we have to have 

w{9) = ah{9) , 
where a is a constant. Next, the bYBe for \Ax{9i)Ai{92)) 

d ff{0)\ 



d9 \gi9) 




whose solution is 



f{9) P-9 



g{9) fi + 9 

for some constant /?. Finally, the bYBe for \Ai{9i)Ai{92)) 



g{9)^lh{9)'- 



9 



(5.43) 

(5.44) 
(5.45) 

|^2(-^i)^2(-^2)) requires 
(5.46) 

(5.47) 

\A2{-9i)Ai{-92)) demands 

(5.48) 
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The reflection matrix is, then, 

m = h{9) ( ^ ) ■ (5-49) 

To obtain the solution where we would have set h{6) = (or equivalcntly w{9) = 0) 
in 5.44, all we should do is to set a = in 5.49. We have checked that 5.49 solves 
all the remaining bYBc's. The function h{6) may be fixed by unitarity and boundary 
crossing-unitarity, but again, we are not going to carry out this computation here. 

One aspect that may look puzzling is the fact that we know that the 0(2) nla model 
is related to the sine-Gordon model, and 5.49 has 3 independent parameters, unlike the 
boundary sine-Gordon model, where one has only 2 independent parameters. The reason 
for that is that we are not looking at the sG model at arbitrary coupling constant, but 
at (3'^ = 8n, and as pointed out in [1], at special points for the coupling constant (like 
this one) there are more solutions than the ones presented there for arbitrary (3. 

We now go back to the solution of unitarity and boundary-crossing symmetry con- 
ditions for the reflection matrix 5.27. The next two subsections refer to this reflection 
matrix. 

5.4 Boundary Unitarity 

The boundary-unitarity condition implies that 

Ro{e)Ro{-e) = 1 , (5.50) 

and 

A{e)A{-e) - B{e)B{-e) = l and A{-9)B{9) + A{9)B{-e) = . (5.51) 

By dividing 5.51 by 1 = R{9)R{—9), and using the explicit forms 5.34, we see that 5.51 
is satisfied trivially. Therefore, all we have to solve from unitarity is 5.50. 

5.5 Boundary Crossing Unitarity 

The boundary-crossing unitarity condition for the diagonal part of the refiection matrix 
is 

?7r ITT IIT 

My - ^) = ((^ - 2)<^i(2^) + ^2(2^) + a3(2e))i?o(y + 9)+ 2ai(2e)A(- + 9) . (5.52) 
Using 5.34 and 5.36, we can rewrite 5.52 as 

RoCy -0)^ RoCy + 0){'Jim + 2a,{29){xC-^ + 9)- 1)) . (5.53) 
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All the other boundary-crossing unitarity equations are equivalent to this one. Using the 
parametreization for c and d introduced after equation 5.36, and after some elementary, 
but tedious, computation, we can cast 5.53 into a quite compact form 



^o(f - e) 



,171 



e 



(5.54) 



Finally we can fix Rq{0) by solving 5.54 and 5.50 



Ro{e) 



U 47r 2-k) V4 "t Att 



^ U 47r 2n)^ U 



+ £)r(| + 



A+2ig 
47r 



X 



2tt^ 



r(| + 



A-2ig 
47r 



2tt>^ \4:' 47r 27r^ 



(5.55) 



where K{6) is the solution of 5.54 and 5.50 without the ,^-dependent factors [15], 



V2 47r 27r/ V-"- 27r/ U 47r 27r/ U 2n ) 



r(i + ^ + f)r(i- 



27r 



M V4 ^ 47r 27r/ ^ U ^ 27r/ 



(5.56) 



This fixes the reflection matrix completely. 

This solution for the reflection matrix reduces to the appropriate diagonal solutions 
in the limits ^ — > (2 Neumann, N — 2 Dirichlet) and ^ — > 00 (all Dirichlet). 

It would be interesting to study the analytical structure of this reflection amplitude, 
but this requires the knowledge of the explicit form of the function ^{g). 



6 Conclusions 

We have found new integrable boundary conditions for the 0{N) nlcr model which are 
nondiagonal and depend on one free parameter. The integrable boundary conditions 
presented here break the 0{N) symmetry to 0(2) x 0{N — 2). A similar conclusion has 
been reached for the SU{N) spin chain in [17]. These boundary conditions, together with 
the diagonal ones proposed in [15], are argued to exhaust the possible integrable boundary 
conditions for the 0{N) nlcr model. We also found solutions of the bYBe which depend 
on one parameter and that can be associated to these boundary conditions, and saw that 
they do satisfy the appropriate limits, which correspond to the coupling constant 51 — > 
or 00. We found the most general solution of the bYBe for the 0(2) nlcr model, which 
is described by a three-parameter family, unlike the well-known two-parameter solution 
of the boundary sine-Gordon. This is understood by noting that the 0(2) nlcr model 
corresponds to the = Stt sine-Gordon model, and as pointed out in[l], at this point, 
and in general for points where X — 8tt/ (5^ — 1, there are additional solutions. 
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In section 4 we established the integrabihty of the nondiagonal boundary conditions, 
by inspecting the spin-4 charges proposed in [9] . There were no constraints on how many 
nondiagonal blocks one can have, whereas we found that there are solutions of the bYBe 
only in the case of one single nondiagonal block. One possible explanation for that is 
that by inspecting higher spin charges there would be further constraints on the possible 
boundary conditions, but this is still to be checked. 

Recently, MacKay and Short [6] have studied the principal chiral model with a bound- 
ary and found an interesting relationship between their boundary conditions and the 
theory of symmetric spaces. Their solutions, though, are quite different from ours, and 
some work should be done in trying to clarify their relationship. 

There are several directions to be pursued now. The S-matrix for the elementary 
excitations in the SO{N) principal chiral and 0{N) Gross- Neveu models, are, up to 
CDD factors, the same as the one for the 0{N) nla model, and therefore the bYBe's 
for all the three models is the same. This means that it should be possible to find 
microscopic lagrangians which provide nondiagonal boundary conditions for these other 
models, similar to the ones presented here, also depending on one free parameter. 

Another problem that arises naturally is to find the exact form of the function ^{g), 
which is, in general quite hard (as for example, in the case of the boundary sine-Gordon 
model). A related problem is to perform a large- expansion for the refiection matrices 
(diagonal and nondiagonal) as a check of their validity. Unfortunately a framework of 
how to perform boundary perturbation theory is yet to be developed. Finally, the study 
of the thermodynamics of the boundary nlcr model would be an interesting, if somewhat 
challenging, problem. 
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